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Q ' Abstract 

C/3 ! 

^ I We introduce a general stochastic model for the spread of rumours, and derive mean- 

• i-H ■ field equations that describe the dynamics of the model on complex social networks 

^. (in particular those mediated by the Internet). We use analytical and numerical 

1-^ I solutions of these equations to examine the threshold behavior and dynamics of 

1 '•^ ' the model on several models of such networks: random graphs, uncorrelated scale- 

free networks and scale-free networks with assortative degree correlations. We show 

►^ i that in both homogeneous networks and random graphs the model exhibits a critical 

00 ' threshold in the rumour spreading rate below which a rumour cannot propagate in 

'^ ■ the system. In the case of scale-free networks, on the other hand, this threshold 

^ , becomes vanishingly small in the limit of infinite system size. We find that the 

I • I initial rate at which a rumour spreads is much higher in scale-free networks than 

f~^ ' in random graphs, and that the rate at which the spreading proceeds on scale- free 

00 , networks is further increased when assortative degree correlations are introduced. 

O ; The impact of degree correlations on the final fraction of nodes that ever hears 

^^ ' a rumour, however, depends on the interplay between network topology and the 

S^ . rumour spreading rate. Our results show that scale-free social networks are prone 

H I to the spreading of rumours, just as they are to the spreading of infections. They are 

■ - - ' relevant to the spreading dynamics of chain emails, viral advertising and large-scale 

information dissemination algorithms on the Internet. 



1 Introduction 



Rumours are an important form of social communications, and their spreading 
plays a significant role in a variety of human affairs. The spread of rumours 
can shape the public opinion in a country [1] , greatly impact financial markets 
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Pl3] and cause panic in a society during wars and epidemics outbreaks. The 
information content of rumours can range from simple gossip to advanced 
propaganda and marketing material. Rumour-like mechanisms form the basis 
for the phenomena of viral marketing, where companies exploit social networks 
of their customers on the Internet in order to promote their products via the so- 
called 'word-of-email' and 'word-of-web' [1]. Finally, rumor-mongering forms 
the basis for an important class of communication protocols, called gossip 
algorithms, which are used for large-scale information dissemination on the 
Internet, and in peer-to-peer file sharing applications 



Rumours can be viewed as an "infection of the mind", and their spreading 
shows an interesting resemblance to that of epidemics. However, unlike epi- 
demic spreading quantitative models and investigations of rumour spreading 
dynamics have been rather limited. An standard model of rumour spreading, 
was introduced many years ago by Daley and Kendall [7115] . The Daley-Kendall 
(DK) model and its variants, such as the Maki-Thompson (MK) model [9], 
have been used extensively in the past for quantitative studies of rumour 
spreading [ T0|ll|12p!3] . In the DK model a closed and homogeneously mixed 
population is subdivided into three groups, those who are ignorant of the ru- 
mour, those who have heard it and actively spread it, and those who have 
heard the rumour but have ceased to spread it. These groups are called igno- 
rants, spreaders and stiflers, respectively. The rumour is propagated through 
the population by pair-wise contacts between spreaders and others in the pop- 
ulation, following the law of mass action. Any spreader involved in a pair-wise 
meeting attempts to 'infect' the other individual with the rumour. In the case 
this other individual is an ignorant, it becomes a spreader. In the other two 
cases, either one or both of those involved in the meeting learn that the ru- 
mour is 'known' and decided not to tell the rumour anymore, thereby turning 
into stiflers [8] . In the Maki-Thompson variant of the above model the rumour 
is spread by directed contacts of the spreaders with others in the population. 
Furthermore, when a spreader contacts another spreader only the initiating 
spreader becomes a stifier. 

An important shortcoming of the above class of models is that they either 
do not take into account the topology of the underlying social interaction 
networks along which rumours spread (by assuming a homogeneously mixing 
population), or use highly simphfied models of the topology pT|IT2] . While 
such simple models may adequately describe the spreading process in small- 
scale social networks, via the word-of-mouth, they become highly inadequate 
when applied to the spreading of rumours in large social interaction networks, 
in particular those which are mediated by the Internet. Such networks, which 
include email networks [Ti|15lll6j . social networking sites [17] and instant mes- 
saging networks [18] typically number in tens of thousands to millions of nodes. 
The topology of such large social networks shows highly complex connectivity 
patterns. In particular, they are often characterized by a highly right-skewed 



degree distribution, implying the presence of a statistically significant number 
of nodes in these networks with a very large number of social connections 



A number of recent studies have shown that introducing the complex topology 
of the social networks along which a rumour spreads can greatly impact the 
dynamics of the above models. Zanette performed simulations of the determin- 
istic Maki-Thompson model on both static [19j and dynamic [20] small-world 
networks. His studies showed that on small-world networks with varying net- 
work randomness the model exhibits a critical transition between a regime 
where the rumour "dies" in a small neighborhood of its origin, and a regime 
where it spreads over a finite fraction of the whole population. Moreno et 
al. studied the stochastic version of the MK model on scale-free networks, 
by means of Monte Carlo simulations [2T] , and numerical solution of a set of 
mean-field equations [22] . These studies revealed a complex interplay between 
the network topology and the rules of the rumour model and highlighted the 
great impact of network heterogeneity on the dynamics of rumour spreading. 
However, the scope of these studies were limited to uncorrelated networks. An 
important characteristic of social networks is the presence of assortative de- 
gree correlations, i.e. the degrees of adjacent vertices is positively correlated 
[23|24fl4|[T5] . Furthermore the mean- field equations used in [22] were postu- 
lated without a derivation. 

In this paper we make several contributions to the study of rumour dynamics 
on complex social networks. First of all, we introduce a new model of rumour 
spreading on complex networks which, in comparison with previous models, 
provides a more realistic description of this process. Our model unifies the 
MK model of rumour spreading with the Susceptible-Infected-Removed (SIR) 
model of epidemics, and has both of these models as it limiting cases. Sec- 
ondly, we describe a formulation of this model on networks in terms of In- 
teracting Markov Chains (IMC) [22], and use this framework to derive, from 
first-principles, mean-field equations for the dynamics of rumour spreading on 
complex networks with arbitrary degree correlations. Finally, we use approx- 
imate analytical and exact numerical solutions of these equations to examine 
both the steady-state and the time-dependent behavior of the model on several 
models of social networks: homogeneous networks, Erdos-Renyi (ER) random 
graphs, uncorrelated scale-free networks and scale-free networks with assorta- 
tive degree correlations. 

We find that, as a function of the rumour spreading rate, our model shows a 
new critical behavior on networks with bounded degree fluctuations, such as 
random graphs, and that this behavior is absent in scale-free networks with 
unbounded fluctuations in node degree distribution. Furthermore, the initial 
spreading rate of a rumour is much higher on scale-free networks as compared 
to random graphs. This spreading rate is further increased when assortative 



degree correlations are introduced. The final fraction of nodes which ever hears 
a rumour (we call this the final size of a rumour), however, depends on an in- 
terplay between the model parameters and degree-degree correlations. Our 
findings are relevant to a number of rumour-like processes taking place on 
complex social networks. These include the spreading of chain emails and In- 
ternet hoaxes, viral advertising and large-scale data dissemination in computer 
and communication networks via the so-called gossip protocols [5]. 

The rest of this paper is organized as follows. In Section 2 we describe our 
rumour model. In section 3 a formulation of the model within the framework 
of Interactive Markov Chains is given, and the corresponding mean-field equa- 
tions are derived. In section 4 we present analytical results for the steady-state 
behavior of our model for both homogeneous and inhomogeneous social net- 
works. This is followed in section 5 by numerical investigations of the steady- 
state and dynamics of the model on several models of social networks: the 
ER random graph, the uncorrelated scale-free networks and the assortatively 
correlated SF networks. We close this paper in section 6 with conclusions. 



2 A general model for rumour dynamics on social networks 



The spread of rumours is a complex socio-psychological process. An adequate 
modeling of this process requires both a correct description of the underlying 
social networks along which rumours spread and a quantitative formulation 
of various behavioural mechanisms that motivate individuals to participate in 
the spread of a rumour. The model described below is an attempt to formalize 
and simplify these behavioral mechanisms in terms of a set of simple but 
plausible rules. 

Our model is defined in the following way. We consider a population consisting 
of N individuals which, with respect to the rumour, are subdivided into ig- 
norants, spreaders and stifiers. Following Maki and Thompson [9j, we assume 
that the rumour spreads by directed contact of the spreaders with others in the 
population. However, these contacts can only take takes place along the links 
of an undirected social interaction network G = (V, E), where V and E denote 
the vertices and the edges of the network, respectively. The contacts between 
the spreaders and the rest of the population are governed by the following set 
of rules 

• Whenever a spreader contacts an ignorant, the ignorant becomes an spreader 
at a rate A. 

• When a spreader contacts another spreader or a stifier the initiating spreader 
becomes a stifier at a rate a. 



In the above, the first rule models the tendency of individuals to accept a ru- 
mour only with a certain probability which, loosely speaking, depends on the 
urgency or credibility of a rumour. The second rule, on the other hand, models 
the tendency of individuals to lose interest in spreading a rumour when they 
learn, through contacts with others, that the rumour has become stale news, 
or is false. In both the Daley-Kendall and the Maki- Thompson rumour mod- 
els, and their variants, stifling is the only mechanism that results in cessation 
of rumour spreading. In reality, however, cessation can occur also purely as 
a result of spreaders forgetting to tell the rumour, or their disinclination to 
spread the rumour anymore. Following a suggestion in [8J, we take this im- 
portant mechanism into account by assuming that individuals may also cease 
spreading a rumour spontaneously (i.e. without the need for any contact) at 
a rate 6. The spreading process starts with one (or more) element (s) becom- 
ing informed of a rumour and terminates when no spreaders are left in the 
population. 



3 Interactive Markov chain mean-field equations 



We can describe the dynamics of the above model on a network within the 
framework of the Interactive Markov Chains (IMC). The IMC was originally 
introduced as a means for modelling social processes involving many interact- 
ing actors (or agents) ^25j. More recently they have been applied to the dy- 
namics of cascading failures in electrical power networks [26] , and the spread 
of malicious software (malware) on the Internet [27]. An IMC consists of A^ 
interacting nodes, with each node having a state that evolves in time according 
to an internal Markov chain. Unlike conventional Markov chains, however, the 
corresponding internal transition probabilities depend not only on the current 
state of a node itself, but also on the states of all the nodes to which it is 
connected. The overall system evolves according to a global Markov Chain 
whose state space dimension is the product of states describing each node. 
When dealing with large networks, the exponential growth in the state space 
renders direct numerical solution of the IMCs extremely difficult, and one has 
to resort to either Monte Carlo simulations or approximate solutions. In the 
case of rumour model, each internal Markov chain can be in one of the three 
states: ignorant, spreader or stifier. For this case we derive below a set of 
coupled rate equations which describe on a mean-field level the dynamics of 
the IMC. We note that a similar mean-field approach may also be applicable 
to other dynamical processes on networks which can be described within the 
IMC framework. 

Consider now a node j which is in the ignorant state at time t. We denote 
with pjj the probability that this node stays in the ignorant state in the time 
interval [t,t + At], and with pj^ = 1 — pjj the probability that it makes a 



transition to the spreader state. It then follows that 

pi = ii- AtA)^ (1) 

where g = g{t) denotes the number of neighbors of node j which are in the 
spreader state at time t. In order to progress, we shall coarse-grain the micro- 
dynamics of the system by classifying nodes in our network according to their 
degree and taking statistical average of the above transition probability over 
degree classes. 

Assuming that a node j has k links, g can be considered as an stochastic 
variable which has the following binomial distribution: 

Il{g,t)=QB{k,ty{l-e{k,t)f-\ (2) 

where 6{k,t) is the probability at time t that an edge emanating from an 
ignorant node with k links points to a spreader node. This quantity can be 
written as 

e{k, t) = Y^ p{k'\k)PMik) ^ E P{k'\k)p'{k', t). (3) 

k' k' 

In this equation P{k'\k) is the degree-degree correlation function, P{sk'\ik) 
the conditional probability that a node with k' links is in the spreader state, 
given that it is connected to an ignorant node with degree k, and p^{k\t) is 
the density of spreader nodes at time t which belong to connectivity class k. In 
the above equation the final approximation is obtained by ignoring dynamic 
correlations between the states of neighboring nodes. 

The transition probability pii{k^ t) averaged over all possible values of g is then 
given by: 



p,,(A;, t) = E (^\ (1 - xMye^k, t)ni - e{k, t)f'^ 

= (l-XAtJ2Pik'\k)p'ik',t)] . (4) 



In a similar fashion we can derive an expression for the probability pss{k,t) 
that a spreader node which has k links stays in this state in the interval 
[t,t + At]. In this case, however, we also need to compute the expected value 
of the number of stifler neighbors of the node at time t. Following steps similar 
to the previous paragraphs we obtain 



Pss{k,t) = [1 - aAtJ2P{k'\k){p'{k',t) + p'{k',t)) 

\ k' y 

X {l-6At). 



(5) 



The corresponding probability for a transition from the spreader to the stifler 
state, psr{k, t) is given by psr{k, t) = 1 — pss{k, t). 

The above transition probabihties can be used to set up a system of cou- 
pled Chapman-Kolmogorov equations for the probability distributions of the 
population of ignorants, spreaders and stiflers within each connectivity class. 
However, ignoring fluctuations around expected values we can also obtain a 
set of deterministic rate equations for the expected values of these quantities. 

Denote with I{k,t), S{k,t),R{k,t) the expected values of the populations of 
nodes belonging to connectivity class k which at time t are in the ignorant, 
spreader or stifler state, respectively. The event that an ignorant node in class 
k will make a transition to the spreader state during [t, t + At] is a Bernoulli 
random variable with probability (1 — pii{k,t)) of success. As a sum of i.i.d 
random Bernoulli variables, the total number of successful transitions in this 
time interval has a binomial distribution, with an expected value given by 
I{k,t){l — pii{k,t)). Hence the rate of change in the expected value of the 
population of ignorant nodes belonging to class k is given by 



I{k,t + At) = I{k,t)-I{k,t) 



1 - ( 1 - \AtY,p"{k',t)P{k'\k)\ 



(6) 



Similarly, we can write the corresponding rate of change in the population of 
spreaders and stiflers as follows 



S{k,t + At) = S{k,t) + I{k,t) 



-S{k,t) 
-6S{k,t) 



1 - ( 1 - XAtJ2p"{k',t)P{k'\k)\ 
1 - ( 1 - aAtY,{p%k',t) + p'{k',t))P{k'\k)\ 



(7) 



R{k,t + At)=R{k,t) 
+ Sik,t) 
+ 6S{k,t) 



1 - ( 1 - aAtY,{p"{k', t) + p''(A;', t))P{k'\k)\ 



In the above equations p^{k,t), p'^{k,t), and p^{k,t) are the fraction of nodes 
belonging to class k which are in the ignorant, spreader and stifler states, 
respectively. These quantities satisfy the normalization condition p*(A;,t) + 
p'^ik, t) + p'^Xk, t) = 1. In the limit At — *> we obtain: 

^^^ = -kV{k,t)Y.p%k',t)P{k'\k) (9) 

k' 



^^^ = kxp\k, t) y: p'{k', t)p{k'\k) 

- kap'{k, t) J2ip'{k', t) + p"{k', t))P{k'\k) - 6p'{k, t). (10) 



^^^ = kap-^{k, t) Y.{p\k', t) + /(A;', t))P{k'\k) + 5p\k, t). (11) 



For future reference we note here that information on the underlying network 
is incorporated in the above equations solely via the degree-degree correlation 
function. Thus in our analytical and numerical studies reported in the next 
section we do not need to generate any actual network. All that is required 
is either an analytical expression for P{k!\k^ or a numerical representation of 
this quantity. 



4 Analysis 

J^.l Homogeneous networks 



In order to understand some basic features of our rumour model we first con- 
sider the case of homogeneous networks, in which degree fluctuations are very 
small and there are no degree correlations. In this case the rumour equations 
become: 



^ = ~Xkp\t)p^{t) (12) 

dp^it) 



\kp\t)p\t) - akp'{t){p\t) + p'it)) 
6p\t) 

akp\t){p\t)+p\t)) + 5p\t), (14) 



dt 

-Sp^it) (13) 

dp^'it) 
dt 



where k denotes the constant degree distribution of the network (or the average 
value for networks in which the probabihty of finding a node with a different 
connectivity decays exponentially fast). 

The above system of equations can be integrated analytically using an stan- 
dard approach. In the infinite time limit, when p''^{oo) = 0, we obtain the 
following transcendal equation for R = p^loo), the final fraction of nodes 
which ever hear the rumour (we call this the final rumour size) 

R=l- e"^ (15) 

where 

' (" + ^)^ (16) 



6 + ak 



Eq. (lisp admits a non-zero solution only if e > 1. For 6^0 this condition is 
fulfilled provided 

^k>l, (17) 

which is precisely the same threshold condition as found in the SIR model 
of epidemic spreading on homogeneous networks [28|29] . On the other hand, 
in the special case 6 = (i.e when the forgetting mechanism is absent) e = 
1 + A/a > 1, and so Eq. (14) always admits a non-zero solution, in agreement 
with the result in [221. 



The above result shows, however, that the presence of a forgetting mechanism 
results in the appearance of a finite threshold in the rumour spreading rate 
below which rumours cannot spread in homogeneous networks. Furthermore, 
the value of the threshold is independent of a (i.e. the stifling mechanism), and 
is the same as that for the SIR model of epidemics on such networks. This 
result can be understood by noting that in the above equations the terms 
corresponding to the stifling mechanism are of second order in p^, while the 
terms corresponding to the forgetting mechanism are only of first order in this 
quantity. Thus in the initial state of the spreading process, where p* ~ and 
p^ ^ 0, the effect of stifling is negligible relative to that of forgetting, and the 
dynamics of the model reduces to that of the SIR model. 



4-2 Inhomogeneous networks 



Next we consider uncorrelated inhomogeneous networks. In such networks the 
degree-degree correlations can be written as [30] : 



P{k'\k) = q{k')-- ^'^^^'^ 



where P{k') is the degree distribution and {k) is the average degree. In this 
case the dynamic of rumour spreading is describe by Eqs. fl MTTj) . Eq. ([9]) can 
be integrated exactly to yield: 



p\k,t) = p\k,Q)e-^^^^'\ (19) 

where p*(A;, 0) is the initial density of ignorant nodes with connectivity /c, and 
we have introduced the auxiliary function 

m = J2q{k) fp\k,t')dt' ^ [\{p\k,t')))dt'. (20) 

In the above equation and hereafter we use the shorthand notation 

{{0{k)))=j:q{k)0{k) (21) 

k 

with 

.(^•) - ^. (22) 

In order to obtain an expression for the final size of the rumour, R, it is more 
convenient to work with 0. Assuming an homogeneous initial distribution of 
ignorants, p\k, 0) = pg, we can obtain a differential equation for this quantity 
by multiplying Eq. (llOp with q{k) and summing over k. This yields after some 
elementary manipulations: 



'^'^ 1 - ((e-^^^*))) -6^-a f'h- {{e-''^^''^))] {{kp\k,t')))dt\ (23) 



dt Jo 

where, without loss of generahty, we have also put pg ~ 1- 
In the limit t ^ oo we have ^ = 0, and Eq. (1231) becomes: 



= 1- ((e-^^<^°°)) -50OO -« / 1 - ((e-^^^(*))) {(kp\k,t')))dt', 

Jo L J 



(24) 



10 



where 0oo = limt_oo0(t). 

For a = Eq. flMl) can be solved explicitly to obtain $00 [28] • For a 7^ we 
solve (^^ to leading order in a. Integrating Eq. (ITUl) to zero order in a we 
obtain 

p«(fc, t) = 1 - e-^'^'^ - 5 r e^(*-*') fl - e-^^'^(*')l dt' + 0{a). (25) 

Close to the critical threshold both 0(t) and (poo are small. Writing (f){t) = 
(poofit), where f{t) is a finite function, and working to leading order in (pac, 
we obtain 

p'{k, t) - -6Xk<p^ t e'^'~''^f{t')dt' + 0{<Pl) + 0{a) (26) 



Inserting this in Eq. fl24l) and expanding the exponential to the relevant order 
in (hn^ we find 



= 0OC [\m) - S - \^{{k^)){l/2 + a((A;))/)0oo] + 0{a^) + 0(0^) (27) 

where J is a finite and positive-defined integral. The non-trivial solution of 
this equation is given by: 



" A2((P))(i + «/((fc)))- 
Noting that {{k)) = {k^)/{k) and ((P)) = {k^)/{k) we obtain: 



2{k){^\-5) 
A2(F)(l + 2a/©)' 



(29) 



(fc) 
This yields a positive value for (p^o provided that 



A (A:) 

Thus, to leading order in a, the critical rumour threshold is independent of the 
stifiing mechanism and is the same as for the SIR model [SSPS] . In particular, 
ioT 6 = 1 the critical rumour spreading threshold is given by Ac = {k)/{k'^), 
and Eq. (!29|) simplifies to: 

^ _ 2(A:)(A-A,) ^g^^ 



^°° A2(A;3)(A, + 2a/)' 
Finally, R is given by 

11 



R = Y,P{k){l-e-''''^^), (32) 



k 



The solution to the above equation depends on the form of P{k). In particular, 
for homogeneous networks where all the moments of the degree distribution 
are bounded, we can expand the exponential in Eq. fl32p to obtain 



^-?^'^'^^^-- A(l^(t;^/) - (^^) 



which shows that i? ~ (A — Ac) in the vicinity of the rumour threshold. For 
heterogeneous networks, one must solve the equation for P{k). This can be 
done for example, as for the SIR model ^28j. 



5 Numerical results 



5. 1 Random graphs and uncorrelated scale-free networks 



We consider first uncorrelated networks, for which the degree-degree correla- 
tions are given by Eq. (ITSl ). We shall consider two classes of such networks. 
The first class is the Erdos-Renyi random networks, which for large N have a 
Poisson degree distribution: 

P{k) = e-^i|^. (34) 

The above degree distribution peaks at an average value (k) and shows small 
fluctuations around this value. The second class we consider are scale-free 
networks which have a power law degree distribution: 

P(A;) = <'^^"' ^"^^^-^ (35) 

otherwise. 

In the above equation fcmin is the minimum degree of the networks and A is 
a normalization constant. Recent studies of social networks on the Internet 
indicates that many of these networks show highly right-skewed degree dis- 
tributions, which could often be modelled by a power-law degree distribution 
[T5|17fl8j . For 2 < 7 < 3 the variance of the above degree distribution becomes 
infinite in the limit of infinite system size while the average degree distribution 
remains finite. We shall consider henceforth SF networks with 7 = 3. 
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Our studies of uncorrelated networks were performed using the above forms of 
P{k) to represent ER and SF networks, respectively. The size of the networks 
considered was A^ = 10^ and A^ = 10^, and the average degree was fixed at 
(k) = 7. For each network considered we generated a sequence of A^ random 
integers distributed according to its degree distribution. The coupled set of 
differential equation fl MTTj) were then solved numerically using an standard 
finite difference scheme, and numerical convergence with respect to the step 
size was checked numerically. In the following and throughout the paper all 
calculations reported are performed by starting the rumour from a randomly 
chosen initial spreader, and averaging the results over 300 runs with different 
initial spreaders. The calculations reported below were performed for networks 
consisting of A^ = 10^ nodes. 

In our first set of calculations we set 5 = 1 and investigate the dynamics as a 
function of the rumour spreading rate A and the stifling rate a. First we focus 
on the impact of network topology on the final size of a rumour, R, which for 
inhomogeneous networks is obtained from 

R = J2p''{k,t^), (36) 

k 

where too denotes a sufficiently long time at which the spreading process has 
reached its steady state (i.e. no spreaders are left in the network). In Fig. 1 R 
corresponding to the ER network is plotted as a function of A, and for several 
different values of the stifling parameter a. It can be seen that in this network 
R exhibits a critical threshold Ac below which a rumour cannot spread in 
the network. Furthermore, just as in the case of homogeneous networks, the 
value of the threshold does not depend on a, and is at Ac = 0.12507. This 
value is in excellent agreement with the analytical results obtained in the 
previous section. We also verifled numerically that, in agreement with our 
analytical flndings, the behaviour of R in the vicinity of the critical point can, 
be described with the form 

i?~A(A-Ac), (37) 

where A = A{a) is a smooth and monotonically decreasing function of a. The 
results are shown in Fig. 2 where R is plotted as function of A for a range of 
values of a, together with the corresponding fits. 

Next we turn to our results for the SF network. In Fig. 3 results for R in 
this network are shown. In this case we also observe the presence of an a- 
independent rumour threshold, albeit for much smaller spreading rates than 
for the ER network. We have verified numerically that in this case the thresh- 
old is approached with zero slope, as can also be gleaned from Fig. 3. Since 
the value of the threshold is independent of a, we can use the well-known 
result for the SIR model (the a = case) to conclude that in the limit of 
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infinite system size the threshold seen in the SF network will approach zero. 
It is therefore not an intrinsic property of rumour spreading on this network. 

In order to further analyze the behavior of R in SF networks, we have numer- 
ically fitted our results to the stretched exponential form, 

i?~exp(-C/A), (38) 

with C depending only weakly on a. This form was found to describe the 
epidemic prevalence in both the SIS and the SIR model of epidemics [3T|28] . 
The results are displayed in Fig. 4, and they clearly indicate that the stretched 
exponential form also nicely describes the behavior of R in our rumour model. 
This result provides further support for our conjecture that the general rumour 
model does not exhibit any threshold behavior on SF networks (at least in the 
limit of infinite systems size) . 

In addition to investigating the impact of network topology on the steady- 
state properties of the model, it is of great interest to understand how the 
time-dependent behavior of the model is affected by topology. In Figs. 5 and 
6 we display, as representative examples, time evolution of, respectively, the 
total fractions of stiflers and spreaders, in both networks for A = 1 and two sets 
of values of the cessation parameters: {5 = 1, a = 0}, and {5 = 0, a = 1}. The 
first set of parameters corresponds to a spreading process in which cessation 
results purely from spontaneous forgetting of a rumour by spreaders, or their 
disinclination to spread the rumour any further. The second set corresponds 
to a scenario where individuals keep spreading the rumour until they become 
stiflers due to their contacts with other spreaders or stiflers in the network. 
As can be seen in Fig. 5, in the first scenario the initial spreading rate of a 
rumour on the SF network is much faster than on the ER network. In fact, 
we find that the time required for the rumour to reach 50% of nodes in the 
ER random graph is nearly twice as long as the corresponding time on the 
SF networks. This large difference in the spreading rate is due to the presence 
of highly connected nodes (social hubs) in the SF network, whose presence 
greatly speeds up the spreading of a rumour. We note that in this scenario 
not only a rumour spreads initially faster on SF networks, but it also reaches 
a higher fraction of nodes at the end of the spreading process. 

It can be seen from Figs, 5 and 6 that in the second spreading scenario (i.e. 
when stifling is the only mechanism for cessation) the initial spreading rate 
on the SF network is, once again, higher than on the ER network. However, 
unlike the previous situation, the ultimate size of the rumour is higher on the 
ER network. This behavior is due to the conflicting roles that hubs play when 
the stifling mechanism is switched on. Initially the presence of hubs speeds 
up the spreading but once they turn into stiflers they also effectively impede 
further spreading of the rumour. 
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5.2 Assortatively correlated scale-free networks 



Recent studies have revealed that social networks display assortative degree 
correlations, implying that highly connected vertices preferably connect to 
vertices which are also highly connected [23]. In order to study the impact of 
such correlations on the dynamics of our model, we make use of the following 
ansatz for the degree-degree correlation function 



P{k'\k) = (1 - P)q{k') + /?4fc'; (0 < /? < 1). (39) 

The above form allows us to study in a controlled way the impact of degree 
correlations on the spreading of rumour. 

Using the above degree-degree correlation function we numerically solved Eqs. 
(iQl fTTj) for a SF network characterized by 7 = 3 and < k >= 7. The network 
size was fixed at A^ = 100, 000, and we used two values for the correlation 
parameter: /3 = 0.2 and j3 = 0.4. Fig. 7 displays i? as a function of A, and for 
a = 0.5, 0.75, 1 (the value of 6 was fixed at 1). 

It can be seen that below A ~ 0.5 a rumour will reach a somewhat smaller 
fraction of nodes on the correlated networks than on the uncorrelated ones. 
However for larger values of A this behavior reverses, and the final size of the 
rumour in assortatively correlated networks shows a higher value than in the 
uncorrelated network. We thus conclude that the qualitative impact of degree 
correlations on the final size of a rumour depends very much on the rumour 
spreading rate. 

Finally, we investigate the effect of assortative correlations on the speed of 
rumour spreading. In Fig. 8 we show our results for the time evolution of 
the total fraction of spreaders, S{t), in scale-free networks consisting of A^ = 
100, 000 nodes and for a correlation strength ranging from f3 = to (3 = 0.4. 
In these calculations the value of A was fixed at 1, and we considered two 
values of a: 0,1. It can be seen that the initial rate at which a rumour spreads 
increases with an increase in the strength of assortative correlations regardless 
of the value of a. However, for a = 1 the rumour also dies out faster when 
such correlations are stronger. 



6 Conclusions 



In this paper we introduced a general model of rumour spreading on complex 
networks. Unlike previous rumour models, our model incorporates two distinct 
mechanisms that cause cessation of a rumour, stifling and forgetting. We used 
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an Interactive Markov Chain formulation of the model to derive determinis- 
tic mean-field equations for the dynamics of the model on complex networks. 
Using these equations, we investigated analytically and numerically the be- 
havior of the model on Erdos-Renyi random graphs and scale-free networks 
with exponent 7 = 3. The critical behavior, the dynamics and the station- 
ary state of our model on these networks are significantly different from the 
results obtained for the dynamics of simpler rumour models on complex net- 
works p!9|20ll21|[22] . In particular, our results show the presence of a critical 
threshold in the rumour spreading rate below which a rumour cannot spread 
in ER networks. The value of this threshold was found to be independent of 
the stifling mechanism, and to be the same as the critical infection rate of 
the SIR epidemic model. Such a threshold is also present in the finite-size SF 
networks we studied, albeit at a much smaller value. However in SF networks 
this threshold is reached with a zero slope and its value becomes vanishingly 
small in the limit of infinite network size. We also found the initial rate of 
spreading of a rumour to be much higher on scale-free networks than on ER 
random graphs. An effect which is caused by the presence of hubs in these 
networks, which efficiently disseminate a rumour once they become informed. 
Our results show that SF networks are prone to the spreading of rumours, just 
as they are to the spreading of infections. 

Finally, we used a local ansatz for the degree-degree correlation function in 
order to numerically investigate the impact of assortative degree correlations 
on the speed of rumour spreading on SF networks. These correlations were 
found to speed up the initial rate of spreading in SF networks. However, their 
impact on the final fraction of nodes which hear a rumour depends very much 
on the rate of rumour spreading. 

In the present work we assumed the underlying network to be static, i.e. a 
time-independent network topology. In reality, however, many social and com- 
munication networks are highly dynamic. An example of such time- dependent 
social networks is Internet chatrooms, where individuals continuously make 
new social contacts and break old ones. Modelling spreading processes on such 
dynamic networks is highly challenging, in particular when the time scale at 
which network topology changes becomes comparable with the time scale of 
the process dynamics. We aim to tackle this highly interesting problem in 
future work. 

M. N. acknowledges the Abdus Salam International Centre for Theoretical 
Physics (ICTP) for a visiting fellowship during which some amendments to this 
work were made. Y. M thanks the support of DGA through a mobility grant 
(MI06/2005) and BT for hospitality. Y. M. is supported by MEC through the 
Ramon y Cajal Program. This work was supported by BT and the Spanish 
DGICYT Projects FIS2004-05073-C04-01 and FIS2005-00337. 



16 



References 



[I] Galam, S. 2003 Modeling rumors: The no plane pentagon french hoax case, 
Physica A 320, 571. 

[2] Kimmel, A. J. 2004 , Rumors and rumor control, J. Behavioral Finance 5, 134. 

[3] Kosfeld, M. 2005 , Rumours and markets, J. Mathematical Economics 41, 646. 

[4] Domingos, P. and Richardson, M. 2002 Mining knowledge-sharing sites for viral 
marketing, in 'Proc. 8th ACM SIGKDD International Conference on Knowledge 
Discovery and Data Mining', FL, USA, pp. 61-70. 

[5] Demers, A., Greene, D., Hauser, C., Irish, W., Larson, J., Shenker, S., Sturgis, H., 
Swinehart, D. and Terry, D. 1987 Epidemic algorithms for replicated database 
maintenance in 'Proc. of the 6th ACM Symposium on Principles of distributed 
computing', Vancouver, Cauda, pp. 1-12. 

[6] Vogels, W., Renesse, R. V. and K. Birman, L. M. 2003 , ACM SIGCOMM 
Computer Communication Review 33, 131-135. 

[7] Daley, D. J., and Kendal, D. G. 1965 Stochastic rumours, J. Inst. Maths Applies 
1, 42. 

[8] Daley, D. J. and Gani, J. 2000 Epidemic Modelling, Cambridge University Press, 
Cambridge, UK. 

[9] Maki, D. P. 1973 Mathematical Models and Applications, With Emphasis on 
Social, Life, and Mangement Sciences, Prentice Hall. 

[10] Lefevre, C. and Picard, P. 1994 , 'Distribution of the final extent of a rumour 
process', J. Appl. Proh. 31, 244. 

[II] Sudbury, A. 1985 , The proportion of population never hearing a rumour, J. 
App. Proh. 22, 443. 

[12] Pittel, B. 1987 On a Daley-Kendal model of rumours, J. App. Proh. 27, 14. 

[13] Noymer, A. 2001 The transmission and persistence of 'urban legends': 
Sociological applications of age-structured epidemic models, J. Mathematical 
Sociology 25, 299. 

[14] Newman, M. E. J., Forest, S. and Balthrop, J. 2002 Email networks and the 
spread of viruses, Phys. Rev. i? 66, 035101. 

[15] Ebel, H., Mielsch, L. and Bornholdt, S. 2002 Scale-free topology of email 
networks, Phys. Rev. E 66, 035103. 

[16] F. Wang, Y.Moreno, and Y. Sun "Structure of peer-to-peer social networks", 
Phys. Rev. E 73, 036123 (2006). 

[17] Csanyi, G. and Szendroi, B. 2004 , Structure of large social networks, Phys. 
Rev. E69, 036131. 



17 



[18] Smith, N. 2002 Instant messaging as a scale- free network, con-mat/0206378 



[19] Zanette, D. H. 2001, Critical behavior of propagation on small-world networks, 
Phys. Rev. E 64, 050901 (R). 

[20] Zanette, D. H. 2001, Dynamics of rumor propagation on small world networks, 
Phys. Rev. E 64, 050901 (R). 

[21] Moreno, Y., Nekovee, M. and Vespignani, A. 2004 , Efficiency and reliability of 
epidemic data dissemination in complex networks Phys. Rev. E 69, 055101R. 

[22] Moreno, Y., Nekovee, M. and Pacheco, A. 2004 , Dynamics of rumor spreading 
in complex networks, Phys. Rev. E 69, 066130. 

[23] Newman, M. E. J. 2003 The structure and function of complex networks, SIAM 
Rev. 45, 167. 

[24] Boccaleti, S., Latora, V., Moreno, Y. Hwang, D.-U., and M. Chavez, 2006 
Complex Networks: Structure and Dynamics, Physics Reports 424, 175. 

[25] Conlisk, J. 1976 Interactive Markov chains, J. Math. Sociology 4, 157. 

[26] Asavathiratam, C, Roy. S, Lesieutre, B. and Verghese, G. 2001 The influence 
model, IEEE Control Magazine 21, 52-64. 

[27] Said, M. R., Oppenheim, A. V. and Lauffenburger, D. A. 2002 , Modelhng 
cellular signal processing using interacting markov chains, in 'Proc. International 
Conference on Acoustic, Speech, Signal processing (ICSB 2002)', Hong Kong. 

[28] Moreno, Y., Pastor-Satorras, R. and Vespignani, A. 2002 Epidemic outbreaks 
in complex heterogeneous networks, Eur. Phys. J. B 26, 521. 

[29] May, R. M. and Lloyd, A. L. 2001 Infection dynamics on scale free networks, 
Phys. Rev. E 64, 066112. 

[30] Pastor-Satorras, R. and Vespignani, A. 2004 , Evolution and Structure of 
the Internet: A Statistical Physics Approach, Cambridge University Press, 
Cambridge, UK. 

[31] Pastor-Satorras, R. and Vespignani, A. 2001 Epidemic spreading in scale-free 
networks Phys. Rev. Lett. 86, 3200. 

[32] Vazquez, A. and Weigt, M. 2003 , 'Computational complexity arising from 
degree correlations in networks', Phys. Rev. E67, 027101. 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



Fig. 1. The final size of the rumor, R is shown as a function of the spreading rate 
A for the ER network of size 10^. The results are shown for several values of the 
stifling parameter a. 




a-K) 



Fig. 2. R is plotted as a function of A — Ac for the ER network of size 10^, using 
different values of a. Solid lines show our numerical fits to the form i? ~ (A — Xc)^, 
with (3 = 1. 
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Fig. 3. The final size of the rumor, R is shown as a function of the spreading rate 
A for the SF network of size 10^. The results are shown for several values of the 
stifling parameter a. 




Fig. 4. R (in log scale) in the SF network of size 10^ is plotted as a function of 1/A 
and several values of a. Solid lines are our numerical fits to the stretched exponential 
form R = B{a) exp(-C(a)/A). 
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Fig. 5. Time evolution of the density of stiflers is shown on the ER (dashed hnes) and 
the SF network (soUd hnes) when the dynamics starts with a single spreader node. 
Results are shown for two sets of values of the cessation parameters {a = 0,5 = 1} 
and {a = 1,6 = 0}. The network sizes are A^ = 10^. 
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Fig. 6. Time evolution of the density of spreaders is shown for the same networks, 
model parameters and initial conditions as in Fig. 5 
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Fig. 7. The final size of the rumor is plotted as a function of A and for several values 
of a in the SF network of size 10^. Results are shown in the absence (solid lines) 
of assortative degree-degree correlations and in the presence of such correlations. 
The correlation strengths used are /3 = 0.2 (short dashed lines) and /5 = 0.4 (long 
dashed lines). 
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Fig. 8. The impact of assortative correlations on time evolution of the density of 
rumour spreaders. Results are shown for the SF networks of size N = 10^ and 
several values of the correlation strength, (5. The upper panel shows results using 
a = and the lower panel those using a = 1. 
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